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I. INTRODUCTION 



The problem of allocating different types of jobs among several 
classes of military manpower is becoming increasingly more important as 
military jobs become more specialized. A model was proposed by Grinold 
[Ref. 1] which constructs a personnel inventory by rank for each of 
several classes of manpower and then allocates that inventory to meet 
billet requirements. The model is designed for long-range planning 
purposes and produces possible inventories based on an optimization 
scheme that sets permitted errors in the allocations. 

The next section reviews the model and offers several points to 
consider concerning its implementation based on practical problems during 
the preparation of this thesis. The review highlights the mathematical 
foundations of the model and defines the essential variables and 
parameters. A complete discussion of the mathematical concepts is found 
in Grinold [Ref. 1] . 

Section III presents a means for measuring the error between require- 
ments and inventories and then expresses this measure as the objective of 
3 types of mathematical programs: linear (goal) programming, elastic 

programming, and quadratic programming. Sections IV, V, and VI present 
a variety of possible formulations for the three programming models from 
which analysis may be accomplished. The formulations presented include 
those proposed by Grinold [Ref. 1] plus several new ones. For the 
quadratic programming models, analytical solutions are derived from which 
computer programs in the APL language were written. 



10 



Section VII is a demonstration of the model using both APL programs 
and a linear programming optimizer available at the Naval Postgraduate 
School. The demonstration is based on typical data from the U.S. Navy 
officer corps. The last section summarizes the results of the investigation 
and proposes areas for further studies. 
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II. REVIEW OF THE MODEL 



A. PURPOSE 

The purpose of the model is to allocate different types of jobs 
among several classes of manpower in a steady state system. The term 
"allocation" refers to a method for distributing jobs among several 
classes of people or people among different types of jobs. The 
allocations are determined by arrays of sharing fractions which are 
defined below. These sharing fractions represent an actual or possible 
management policy. 

The term "steady state" refers to the equilibrium condition prevailing 
in the model. This condition arises from the assumption of a fixed 
longitudinal manpower flow in each manpower class used for the analysis. 
This concept is defined for the model in Grinold [Ref. 1], and is discussed 
in detail in Grinold and Marshall [Ref. 2], and in Bartholomew and Forbes 
[Ref. 3]. Further discussion of the steady state conditions is found in 
the section which defines the steady state variables. 

The final item of interest in stating the purpose deals with the 
planning aspects of the model. Since this analysis is done for steady 
state conditions, it can be used most effectively to test current or 
alternate management policies for their long-term effect and provide 
answers to the "what if" policy questions. As an aggregate planning 
tool, the model is intended to test long-term policies that may be used 
later for daily operation of a manpower system. Since conditions of 
steady state may be significantly different from current conditions. 
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the analyst should be aware of the limitations of the model. Often, 
the relative differences between competing policies or the magnitude of 
the effect of a policy may be the point of interest for which an analysis 
is being performed. The model should be used with these considerations 
in mind. 

The model is motivated by a study of the U.S. Navy officer corps 
and examples from this manpower system will be used when explaining 
concepts as an aid to better understanding of the model. A description 
of the structure of the U.S. Navy officer corps for the unrestricted line 
officers is found in NAVPERS 15197A [Ref. 4]. 

B. DEFINITIONS AND TERMINOLOGY 

1. Attributes 

The definition and understanding of variables in the model is 
complicated by the large number of variables present. The most important 
variables relating to concepts are presented in this section. Other 
variables relating to specific optimization techniques or used for 
mathematical clarity will be introduced as the need for them arises. A 
list of variables used in the model is found in Appendix A. 

The requirements and the inventories possess three common attributes. 
They are: 

• MANPOWER CLASSES--the types of manpower; indexed by k=l,2,...,K. 

• CAREER STA6ES--the experience levels required; indexed by 

i=l,2,...,I. 

• JOB CLASSES--the types of jobs to be filled; indexed by j=l,2,...,J. 

The clear definition of attributes is closely related to the purpose 
of analysis. In the demonstration which follows, the manpower classes 
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represent types of officers (e.g., Pilot, Naval Flight Officer); the 
career stages represent ranks (e.g., Lieutenant, Commander); and the job 
classes are naval officer billet categories (e.g., 1310, 1320). 

The level of detail is also closely related to the purpose 
of analysis. For example, the manpower class of pilot could be further 
broken down into pilots by type of aircraft, but this imposes increased 
data requirements which may not be needed if the purpose of analysis is 
to study Navy-wide trends rather than issues specific to the Navy's 
aviation community. 

The attributes may apply to either requirements or inventories. 

In Navy terminology, requirements are usually stated in terms of "billets." 
Some confusion may arise since this is synonymous with the term "job." 

The reader should bear in mind that requirements may be referred to as 
billets or jobs which are distinct from the attribute of "job type." 
Likewise, the term "inventories" is used interchangeably with stocks, 
personnel, and people, but is distinctly different from the term "manpower 
types." The use of these terms should be clear from the context in which 
they occur. 

The demonstration highlights another problem found in practical 
application. The types of people and types of jobs are often expressed 
by the same terminology. For example, 1310 is both an officer designator 
and a billet designator for the manpower class "pilot" and for the job 
class "pilot." To avoid confusion, the word equivalent "pilot" will be 
used when speaking of 1310 officers and the billet designator "1310" will 
be used when speaking of a corresponding type of job. 
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other interpretations of the model based on current military 
personnel problems are possible. For example, manpower classes could be 
defined as enlisted personnel by mental group category and the purpose of 
analysis could be to study the effects on billets of increasing require- 
ments for the higher mental groups with a shrinking manpower pool from 
which to provide accessions. Another possibility would be to define 
manpower classes by sex and job classes by types of sea or shore billets 
and then study the implications of increased accession of women on 
sea-shore rotation cycles. 

2. Steady State Variables 

Let n. be the number of time periods an individual must remain 
in the organization before entering stage (i+1). For example, naval 
officers in the rank of lieutenant (jg) will remain in the Navy four 
years before being promoted to the rank of lieutenant, thus n^. for 
lieutenant (jg) is four. The total amount of time a person may spend in 
stage i is: 

"i - "i-l ("o ' 

An example of the amount of time spent in each career stage of the Navy 
officer corps will be given in Section VI (Fig. 4) for the demonstration. 
The steady state variables are defined as follows: 

t Survivor fractions by rank and manpower class--this is the matrix 
S indexed over (I,K). The (i,k)^*^ element, s.j^, represents the 
fraction of type k people who remain in the organization more than 
i stages. The element s-j^ may also be interpreted as the 
probaoility that a type k person remains in the organization 
longer than n^ time periods. 
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• Expected stage length is the matrix W whose (i,k)^ element 
represents the amount of time a person in manpower class k is 
expected to spend in stage i. W may also be defined as the 
product of the probability of a type k person reaching stage i and 
the length of stage i: 

"ik ' 

A more complete discussion of survivor fractions and expected stage 
lengths may be found in Grinold [Ref. 1] and Grinold and Marshall [Ref. 2]. 

3 . Requirements 

Requirements are expressed in three ways: 

• Billets by rank and job type--the number of jobs at each stage to 
be filled. This is a matrix B whose elements are indexed by 

• Billets by rank and people type--the number of people required at 
each stage. This is not the normal way in which requirements are 
stated and should not be confused with personnel inventories. 

This is a matrix P whose elements are indexed by (i,k). 

• Target al locations--the number of (i,k) people desired in (i,j) 
jobs. This is a three-dimensional array T with elements indexed 
by (k,i,j). 

4. Inventories and Accessions 

Finding optimal inventories and accessions will normally be the 
goal of this model. They are defined below. 

t h 

• Accessions are denoted by a vector Y whose k ' element represents 
the number of people entering manpower class k each year. 
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Inventories are expressed in three ways: 

• Inventories by rank and job type--the matrix X whose elements are 
indexed by (i,j) defines the number of jobs at each stage actually 
filled. This is not the normal way in which inventories are 
stated and should not be confused with billet requirements. 

• Inventories by rank and people type--the actual number of available 
people by manpower class and rank. This is the matrix Z whose 
elements are indexed by (i,k). 

• Actual al location--the actual number of (i,k) people who will fill 
(i>j) jobs. This is the three-dimensional array A with elements 
indexed by (k,i,j). The notation for requirements and inventories 
is summarized in Table 1. 



TABLE 1 

REQUIREMENTS AND INVENTORIES 





Job Type 
and Rank 


People Type 
and Rank 


People Type, Rank 
and Job Type 


Requirements 




Pik 


\ij 


Inventories 


''ij 


"ik 


\ij 



5 . Sharing Rules 

There are two rules used to distribute the (i,j) jobs among the 
(i,k) manpower types or vice versa. These rules are known as job sharing 
and people sharing. 
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f Job sharing is defined by a three-dimensional array F whose 

element gives the fraction of (i,j) jobs that should 
be performed by manpower type k. 

• People sharing is defined by a three-dimensional array G whose 
(k,i,j) element gives the fraction of (i,k) manpower types who 
should perform job type j. 



C. DETERMINATION OF REQUIREMENTS 

The target allocation (T), which is the number of type j jobs to be 
filled by class k people of rank i, is determined in one of two ways: 

1, By specifying the desired billets by people type and rank (P) and 
the rule for sharing these people among the different billets (G). 
Note that G times P distributes the requirement for the k^ class 
of people in rank i over the job class j and produces the target 
allocation T: 

^kij Pik " '^kij' 

2. By specifying the billets by rank and job type (B) and the rule 
for sharing these billets among the different types of people (F). 
Note that F times B distributes the j^^ type of job in rank i 
over the people class k and also produces the target allocation 

T: 



^kij ^ij " ^kij* 

In order to maintain consistency, the following equation must hold: 

^kij Pik ■ ^kij " ^kij ^ij ^ 

The argument from which this type of equation is derived will be 
referred to as the "consistency rule." Furthermore, all people classes 
must be filling some job and all job classes must be filled by some 
people which means that the equation: 
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implies that 



z 

j 



^kij 



= 1 



Pik 



2 f, . . b . . 

j 'J 



and the equation: 



J 'kij = 1 



implies 



b . . = z q, . . p 
ij 1 ^ ^kij ^ik 



( 2 ) 



( 3 ) 



Thus, one may either start with F and B and calculate P from equation (2) 
and then G using equation (1), or start with G and P and calculate B from 
equation (3) and then F from equation (1). 

For the Navy example which follows, it should be noted that the usual 
manner of determining target allocations will be to start with a given 
requirement of job types by each rank (B) derived from the force structure. 
Since it is more common to think of people being assigned to jobs rather 
than jobs to people, a typical sharing policy is more easily understood 
in terms of distributing people over the required job types (G). Thus, 
it would be helpful to be able to calculate F and P based on B and G. 
Unfortunately, this cannot be done directly. One solution to this 
problem is a least squares approach using equation (3). First determine 
the least squares fit values for P, and then solve for F using equation 
(1). The least squares objective is: 



mi n 
P 



z z 



i j ^k 



(: ^kij Pik 




2 . 
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Notice that this minimization is accomplished by minimizing each of 
the following I expressions: 



For a particular i, the above can be written in matrix notation; 

min (G'P - B)' (G'P - B). 

Thus, the i^ column of P is given by the optimal solution: 

P = [G G'] G B . 

Now calculate f using equation (1) to obtain: 

-1 



^kij " < 



^kij ^ik 
0 



if b.j / 0, 
if b.. -0. 



D. DETERMINATION OF INVENTORIES 

The actual allocation (A), which is the number of people of type k in 
rank i who will be available to fill job type j, is determined using the 
same consistency rule invoked for equation (1); 

^kij ^ik " ^kij " ^kij ^ij 

and the implications similar to equations (2) 
hold. The equation: 



(4) 



and (3) above must also 



implies that 



and the equation: 



t g, . . = 1 

j "’J 



"ik = ' 'kij ' 

vJ 



I hu - 1 



impl i es 






(5) 



( 6 ) 
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In addition, the model is constrained by the steady state requirement 
that the expected inventory of type k people in rank i is equal to the 
number of accessions of type k people times the expected duration of type 
k people in rank i. In equation form, this constraint is: 



Since inventories may be computed from accessions (Y) and expected 
stage lengths (W), the optimization techniques to determine inventories 
found in the next sections will usually have Y as the decision variable. 
If the inventories (X or Z) are the decision variables, then a method is 
needed to determine the corresponding accessions. From equation (7), 
note that for each there are I equations. Once again, a least 
squares approach may be used to find the best fit values of yj^ for this 
system of equations. The objective is: 



min 

y 



(w,- 



ik 



■ 'ik>‘ 



This yields the result: 



^k = '"ik ^-k) for k=l,2,...,K. 
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III. OPTIMIZATION TECHNIQUES 



A. PURPOSE 

This section discusses the use of mathematical programming to deter- 
mine desired values for the variables defined in Section II, The degree 
to which these techniques may be used depends on availability of data, 
level of detail in attribute definition, and desired accuracy of results. 
The scope of the problem to be solved determines the type of the optimi- 
zation technique and the way in which the problem is formulated. 

The purpose of this section is to present a method by which the 
quality of an allocation may be measured. Such a measure will be used as 
the basis for the objective functions for three types of mathematical 
programs: linear (goal), elastic, and quadratic programming. Formula- 

tions for each of these will be presented in subsequent sections. 

B. FORMULATION OF GOALS 

In formulating the optimization, special care must be taken in 
determining the objective to be achieved and the decision variable to be 
determined. The overall objective is to match actual allocations to 
target allocations with respect to a decision variable. The methods 
described here seek to do this by minimizing the weighted percent differ- 
ence between inventories and requirements. The decision variable to be 
determined may be any one of the previously defined variables. All other 
variables are then treated as parameters. 

Listed below are the possible decision variables and a brief discussion 
of how they may affect the optimization. 
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t Requireinents--Errors may be minimized with either billets by job 
type and rank (B) or billets by people type and rank (P) as the 
decision variable. Normally B is treated as a parameter derived 
from the force structure. For example, if ten aviation squadrons 
of a certain type are required by the force structure, then 
various numbers of officer pilots of different ranks will be 
needed to fill the aviation type billets found in those squadrons. 
On occasion, it may be of interest to allow the requirements to be 
decision variables and check the resulting optimal values given 
fixed inventories and sharing rules. 

0 Inventories and Accessions--The decision variable may be either 
the inventories or accessions. Due to the steady state assumption, 
it is sufficient to know the accessions and the expected length of 
time in each career stage in order to calculate expected inven- 
tories. For this reason, the accessions (Y) will normally be the 
decision variables. This approach requires that billets (B or P), 
sharing rules (F or G), and expected duration in each stage (W), 
be treated as parameters. The objective functions in Grinold 
[Ref. 1 ] are of this type. 

0 Expected Duration in Each Career Stage--This variable may be 

the decision variable if billets, accessions and sharing rules are 
known. If certain assumptions are made about the conditions from 
which an optimal W is derived, then it is possible to relate the 
optimal W to a survivor function. 

0 Sharing Rules--Minimi zing errors with F or G as the decision 
variable will result in the optimal "mixing" policy. If the 
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objective is to find the optimal policy based only on billets, 
inventories and expected stage lengths, then this approach is 
reasonable. 

C. MEASUREMENT OF ERROR 

As stated previously, optimization is performed to minimize the 
weighted error between requirements and inventories. The error for 
the linear models and quadratic models may be expressed as the excess or 
deficit of the three quantities shown in Table 2. 



TABLE 2 

TYPES OF ERROR 



Itkii - : the error between targets and 

^ ^ allocations by job type, rank, 

and people type. 

■ Pii^l • the error between the inventory 
of (i,k) people and the corres- 
ponding billets. 

e- • = |b. . - X . . I : the error between the (i,j) billets 

and the corresponding inventory. 



The first error shown involves all three attributes and may be 
obtained from either the second or the third one since these errors may 
be distributed using the job sharing or people sharing rule: 

^kij ‘^ik " l^kij ’ \ij' " ^kij ®ij 
where d^^ and e^j are defined in Table 2. 

For this reason, the last two types of errors, those by rank and job 
type and those by rank and people type, will be used in the following 
discussion. In the linear programming model, it is important to know if 
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the error is a shortfall or overfill. In keeping with the notation in 
Lee [Ref. 5], shortfalls will be denoted by the superscript minus and 
overfills by the superscript plus. 

Weighting factors are used to allow trade-offs between billets of 
differing attributes. The weighting factors are based on total number of 
billets to be filled and the relative importance of billets. A shortfall 
of twenty billets of a particular type will have a different impact 
depending on whether there were 100 or 1,000 billets to be filled. The 
error values will have a common metric if they are converted to percent 
errors. Thus, a better measure of the error would be (z-p)/p or (x-b)/b. 
However, the actual measure may still differ depending on the nature of 
the job type, rank, or personnel type. For example, a five percent 
shortfall of pilots who are of the rank captain may not be as serious as 
a five percent shortfall in pilots who are lieutenants. Thus the percent 
shortfall in meeting a critical target should be weighted more heavily 
than that of a less critical target. In order to handle this problem, 
additional parameters indicating the percent trade-off between ranks, job 
types, and people types must be introduced to the model. One way in 
which this may be done is to specify parameters which will serve as 
benchmarks to measure the permitted range of unit error in overfilling or 
underfilling a requirement. 

For example, let a five percent unit error be the permitted shortfall 
of 1310 (pilot) lieutenant jobs. Now for another job, say 1310 captains, 
the relative trade-off for an underfill may be found by answering the 
question, "What percent under target in the assignment of 1310 captain 
billets is as serious as a five percent shortfall in filling 1310 
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lieutenant billets?" In this way, essential judgments about the trade-off 
between differing ranks and job types may be made. The same process 
could also be done for ranks and people types or the consistency equation 
could be used to produce the trade-off parameters. The notation for these 
parameters is shown in Table 3. These trade-offs may now be added to the 
penalty function by multiplying the percent error by the inverse of the 
permitted unit error. Thus the actual weighting factors for each error 
is the inverse of the permitted unit percent error times the requirements. 
Let u and v be the weighting factors for job types and people types, 
respectively. As before, let the superscript plus denote weights associ- 
ated with overfills and the superscript minus denote weights associated 
with shortfalls. The notation for these parameters is in Table 4. The 
total penalty value is now defined as the weighting factor times the 
error. This penalty value is a measure of error that takes into account 
both the varying sizes of the different requirements and their relative 
importance. The notation for the penalty measure derived thus far is in 
Table 5. 

Obtaining valid weights is of particular importance to the model 
since the final result of the optimization could be very sensitive to 
these parameters. There are several ways in which they can be estimated. 
Expert judges could arbitrarily select either the parameters to be used 
or the range of minimal values for excess and deficit. For example, the 
requirement for a particular job type and rank may be 100. If a permitted 
range of error is 95 to 110, then (> = 5 percent and 9 = 10 percent. The 
range of error is a useful concept and the range bounds F ^ and £ are 
defined in Table 6. If the lower unit errors are greater than 100 
percent, then it would be possible to have a negative value for the lower 
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TABLE 3 






PERMITTED UNIT PERCENT 


ERROR 




by job type 


by people type 




and rank 


and rank 


shortf a1 1 


ij 


^•k 


overf i 1 1 


9. . 
TJ 


^k 



TABLE 4 

WEIGHTING FACTORS 





by job type 


by people type 




and rank 


and rank 


shortf al 1 


U-. = {'!»•• b--) ^ 
ij ^ ij 


'^ik ~ ^'^ik Pik^ 


overf i 1 1 


u"!". = (9. . b. 
ij ^ iJ 


'^ik " ^*ik Pik^ 



TABLE 5 

WEIGHTED PENALTY MEASURE 

by job type by people type 

and rank and rank 



shortf a1 1 



overf i 1 1 



iJ 





- _ Pik ■ ^ik 

V . , d . , = 

Ik 'k 5,^ p,|^ 



>0 ^0 



+ ,+ 
''ik ^-k 



^ik ' Pjk 
*ik Pik 
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TABLE 6 



ERROR RANGE OF REQUIREMENTS 







by job type 
and rank 






by people type 
and rank 


lower 


error bound 


(1 - tij) b,j. 




(1 


- *ik> Pik ' £,'k 


upper 


error bound 


(1 * 9,j) 




(1 


♦ik* Pfk ■ ^'k 



bound. To avoid this possibility, it will be required by definition that 
the lower bounds be greater than or equal to zero. It is also possible to 
estimate the permitted unit error by job type if given people type and 
vice versa using the consistency rule resulting in the equations: 



— ik ^kij 




fkij 


all 


k, i,j; 


^ik ^kij 


' ^ij 


^kij fp-- 


all 


k, i,j. 



D. OBJECTIVE FUNCTIONS 
1 . General 

Let E.. be the penalty value function of errors in job class j 

* J 

and rank i, and let D.j^ be the penalty value function of errors in 
people class k and rank i. Then three possible mathematical programming 
objective functions would be: 

(1) min z z E-., i.e., minimizing error by job types; 

i j ^ 

(2) min z z D-. , i.e., minimizing error by people types; 

i k 

(3) min z z E.. + z z D., , i.e., minimizing error by job types 

i j i k 

and people types. 

Minimizing the third objective function should produce a close match between 
targets and actual allocations and was used in the mathematical programs in 



28 



Grinold [Ref. IJ. The use of objective (3) requires more data input 
for parameters than either objectives (1) or (2). Objective functions 
based on each of the three representations will be used in the following 
sections . 

In order to shorten mathematical notation in the following sections 
on specific objective functions, it will be assumed that the decision 
variable of interest is the inventory X and the objective function is 
type (1) from above. Also, the subscripts will be dropped for clarity. 

2 . Linear (Goal) Programming 

One objective function that measures the weighted percent errors 
is a piecewise linear convex function that has a value of zero if the 
actual allocation meets the target, and increases at a rate inversely 
proportional to the permitted error in meeting that target. This means 
the function is unity when the decision variable assumes either the lower 
or the upper error bound. Therefore, the function will satisfy the 
following three conditions: 



(i) 


E(x) = 0 


if X = b , 


(ii) 


E(x) = 1 


if X = ^ , 


( i i i ) 


E{x) = 1 


if X = 6 . 



Such a function, as shown in Figure 1, is: 



E(x)=™ax 

= max (u'*’ e'*’, u" e”) . 

Notice that the penalty value will be relatively small (less than 1) for 
values between ^ and b and that it will increase at a rate l/(b '(') if 
X < b or l/(b 9) if b < x. 
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E(x) 



1 




b 



b F 



Fig. 1. Piecewise Function for Linear Programming, 



3 , Elastic Programming 

For this type of programming model, the definition of the error 
will be modified slightly as shown in Table 7. The weighting factors 
will be changed since the actual error is now slightly different due to 
the use of the error bounds rather than the actual requirements . In 
order to rotational ly indicate this difference, the use of underscore and 
overscore for weights and errors will be used in place of superscript 
plus and minus. This type of programning can use a piecewise linear 
objective function of the following type: 



This function is shown in Figure 2. Notice that the penalty is 0 for 
values of x within the range of permitted errors (^, b) and increases at 
a rate l/(^ '!») if x is less than ^ or 1/(F 9) if x is greater than F. Also, 
notice that if ^ = b = b then this objective is the same as the goal 






- x) 



if X < b , 



E(x) 




'!» b 
0 



if b < X < b , 



if b < X , 



9 F 

= max (u £, 0, F F) . 
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TABLE 7 



ERROR FOR ELASTIC PROGRAMMING 

by job type by people type 

and rank and rank 



shortfall 




lik 


” -2.ik 


- ^ik 


overf i 1 1 


■ Xij - 


^ik 


II 

tsi 


- Pik 



E(x) 

X 

k 6 

Fig. 2. Piecewise Function for Elastic Programming 

programming objective. Thus, the goal program is a special case of the 
elastic program. 

4. Quadratic Programming 

Several quadratic programming formulations are presented in 
Section VI. These formulations are easy to use and give quick results 
for formulations that can be expressed without inequality constraints. 

The quadratic objective function to be used must have the same properties 
mentioned for the piecewise linear function and it must also be greater 
than zero. These properties may be stated mathematically as: 
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(i) 


E(x) 


= 0 


if X = b , 


(ii) 


E{x) 


= 1 


if X = ^ , 


(iii) 


E(x) 


= 1 


if X = F . 


(iv) 


E{x) 


> 0 





Unfortunately, quadratic functions which use the previously defined 
penalty measure must have symmetric penalties in order to satisfy the four 
conditions. That is, < 1 ) must equal 9 and 5 must equal 4 > when considering 
the decision variable z. The following objective satisfies these 
conditions : 

E{x) = ^ when e = . 

This type of penalty function is shown in Figure 3. An approximation 
that is useful if e and are close but not equal is described in 
Grinold [Ref. 1]. 



E{x) 

I 

I 

I 

I 

I 

1 1 



^ b b 

Fig. 3. Function for Quadratic Programming. 



E. PROGRAMMING NOMENCLATURE 

It will be useful to have a notational scheme to identify the 
various programming types, decision variables, and objective functions 
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found in the following sections. Let LP, EP, and QP refer to linear 
(goal), elastic, and quadratic progr arming, respectively. Let the 
objective function that sums errors by job class be a type 1 objective 
function, the objective function that sums the errors by people class be 
type 2, and the one that sums errors by both job class and people class 
be type 3. Thus, the mathematical program LPl(Y) is a linear program 
that uses the first objective function with Y as a decision variable. 
Similarly, QP2(Y,X) is a quadratic program with the second type of 
objective function that uses Y and X as the decision variables. 

F. ADDING COST CONSIDERATIONS 

Grinold [Ref. 1] also discusses a way in which cost considerations 
may be added to the model. Let C|^ be the annualized cost of a type k 
accession. Then the multi-attribute objective function: 



where 0 _< x £ 1, will minimize the objective with respect to both cost 
and weighted penalties. This objective leads to a family of objective 
functions depending on the value of the parameter x. Notice, however, 
that the units of the objective function are now in terms of cost units 
and percent error units. 

The use of such multi-attribute functions is advantageous when 
the annualized cost of one type of manpower is significantly different 
than another type. For example, the annualized cost of a general unre- 
stricted line officer is vastly less than that of a pilot due to training 
costs, pay bonuses, etc. 
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IV. LINEAR (GOAL) PROGRAMMING FORMULATION 



A. GENERAL 

A linear programming formulation for minimization of errors is 
found in Grinold [Ref. 1] and is derived from the goal programming 
concept as described in Lee [Ref. 5] and Hillier and Lieberman [Ref. 6]. 
Goal programming is a modification and extension of linear programming 
that allows a simultaneous solution of a system of complex objectives 
rather than a single objective. Notice that the goal of minimizing 
errors by differing ranks, job types, and people types constitutes a 
system of conflicting objectives. A solution which forces exact achieve- 
ment of goals for the lower ranks may result in extreme deficits in the 
higher ranks while one which exactly meets requirements in the higher 
ranks may result in large excesses in the lower ranks. Thus, the goal 
approach is used to seek the best possible solution given the stated 
trade-offs. 

In implementation of the following programs, the user should be aware 
of the tactical problems in solving LP's. It may be advantageous in 
some cases to solve the dual formulations of the problems depending on 
the number of variations and the number of constraints. An example of 
forming the dual for LP3(Y) and LP3(X,Y) is in Grinold [Ref. 1] and a 
general discussion of the topic is in Luenberger [Ref. 7]. 

B. ACCESSION AND INVENTORIES 

The first formulation will be designated LPl(Y) ana is shown in 
Figure 4. Recall that LPl formulations involve errors of rank and 
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min 






Efu. 

jLij 



e+ + u“ s” 1 
ij ij ijJ 







g w y ♦ 
kij ik k 




e* 

ij 



y 




e-*" > 0 

ij 



= b for all i,j 
ij 

for all k,i,j 



Fig. 4. LP1{Y) 

job type. The constraint arises from the conservation of flow derived 
from equation (6) and sharing property in equation (7). An alternate 
formulation that considers both types of errors is LP3(Y) as shown in 
Figure 5. This is labelled in Grinold [Ref. 1] as LP-P. 
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Figure 6 displays the formulation which uses both X and Y as the decision 



lies. It is 


LP3(X,Y) 


and 


is found in 


Grinold [Ref, 


•-1] ' 


as LP-B. 


lin ^ ^ 


lu* 9 + 
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a- e* 


> 


E Erv*" d+ ♦ 


r- d-1 


i j 


Lij ij 




ij i j 
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i k L ik 


ik 


ik ikJ 




s. t. z 
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j 


9 “ 

ij 


ij 


= b for 

ij 


all 


i» j 




w y 




d- - 


d+ 


= p for 


all 


if k 




ik 


k 


ik 


ik 


ik 








E 


f 


z 


- w 


7 =0 


for 


all i,j 




j 


kij ij 




ik k 






^3' 


k 1 


d' 


9t ./ 

13 ik 


d-** > 0 for all k,i,j 

rk 



Fig. 6. LP3(X,Y) 



C. EXPECTED STAGE LENGTHS 

Another application using goal programming would be to minimize 
the errors with respect to expected stage lengths (W). Recall that the 
stage length for rank i, person type k is the expected amount of time a 
recruit in class k will spend in rank i. An added constraint for this 
formulation is that the stage lengths have upper and lower bounds. In 
the Navy officer example, such bounds would reflect navy policy or legal 
restrictions imposed by Congress. Let these upper and lower bounds be 
defined by matrices ^ and W, respectively. LPl(W) may be expressed as 
shown in Figure 7. 
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min 



^ Z-lu+ 6 + + U“ 6 " 1 
i jLij ij ij ijj 



s. t. 



E 

Ic 






€- , €+ > 0 
13 ij 



i < H < U 



s b for all i,j 
ij ij 

for all lc,i,j 



Fig. 7. LPl(W). 



D. SHARING RULES 

Of particular interest is a mathematical program to specify the 
sharing fraction that represents the policy for distributing job classes 
among the people types (F) or people classes among job types (G). An 
additional constraint on this formulation is that F must sum to one over 
the index k and G must sum to one over the index j. Also, the sharing 
fractions may be bounded above and below. Using the same convention as 
used in the case of stage lengths, let IF and G be the lower bounds and 
let T and ^ be the upper bounds on F and G, respectively. Notice that if 
no sharing is possible, specifying that both the upper and lower bound be 
zero sets the value at zero. LPl(G) and LP2(F) are formulated as shown 
in Figures 8 and 9. 
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j 



e- , e+ >0 for all i,j 

ij ij 

G < G < G 



Fig. 8, LPl(G) 




Fig. 9. L?2(F) 
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E. INTERPRETATION FROM LP THEORY 



In classical LP terms, the errors e^, e", d^, and d" represent logical 
variables introduced to maintain the equality constraints. The errors e~ 
and d’ may be thought of as artificial variables and the errors e^ and d^ 
as surplus variables and hence, they measure the distance from a particular 
solution to the requirement. These objective functions are mathematically 
equivalent to a Lagrangian form with cost coefficients of the decision 
variables equal to zero (see Duff [Ref. 8: p. 64]). The weighting factors 
u^, u’, v^, and v~ represent bounds on the variables of the dual formulation. 
In fact, an interpretation of a particular weighting factor is that it 
defines the penalty cost per unit violation of its associated constraint. 
Therefore, weighting factors which are approximately equal have associated 
constraints which are of the same relative importance, if binding. In 
other words, the incremental penalty for violating each such constraint 
is approximately the same. Such information should be of use in analyzing 
the results of these programs. 
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V. ELASTIC PROGRAMMING FORMULATION 



A. GENERAL 

The elastic programming model is described in Duff [Ref. 8]. The 
formulation is quite similar to that of the goal program in the last 
section. However, there is now a subtle but important difference. In 
the goal program, the objective was to allocate people as closely as 
possible to an exact number of billets (either B or P). In the elastic 
program, the objective is now to allocate people so as to be within a 
specified range of the billets to be filled. For this reason, the actual 
error will be indicated by a change in the notation to the use of under- 
score and overscore in place of superscripts plus and minus. 

Before starting the formulations, it is necessary to define a variable 
to represent the slack from the upper bounds (B or P). Let q. . and r.. 

I J 1 K 

represent these slack variables for the (i,j) and (i,k) inventories. Also 
recall that £, u, v are the weighting factors and e, ¥, d” are the 
errors from the desired ranges. 

B. ACCESSION AND INVENTORIES 

The first formulation shown in Figure 10 below is similar to the 
first formulation in the previous chapter and is denoted EPl(Y). The 
variable Y once again represents an upper bound on the accession vector Y. 
An alternate formulation that considers both type of errors is EP3(Y) 
presented in Figure 11. 
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Fig, 10. EPl(Y) 



min 
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kij ik k 



+ u jd +? ^1 
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♦e -¥ +q =F for all i,j 
“ij ij ij ij 



w y+d - T +r =p for all i,k 
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0 < Y < y 



0 < c i B - B 
0 < fl < ? - P 

3 /i/® »s -0 for all k,i,j 

ik ik ij ij 



Fig. 11. EP3(Y) 



C. EXPECTED STAGE LENGTHS AND SHARING RULES 



Elastic programming may also be used to determine expected stage 
lengths or sharing rules. Only one example, EP2(F) in Figure 12, will be 
presented since the formulations are quite similar to those in the last 
section. Recall that an 'additional constraint exists for the sharing 
rules since they must sum to one over the appropriate index. An example 
of this is EP2(F) shown below. 




D. INTERPRETATION FROM LP THEORY 

As was the case for the linear (goal) program, the elastic program 
is mathematically equivalent to a Lagrangian form. The variables £ and ^ 
are artificial and measure the distance from a particular solution to the 
lower error bound. The surplus variaoles "e and measure the distance 
above the upper error bound and the slack variables q and r measure the 
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distance from a solution to the associated upper bound. The weighting 
factors continue to represent bounds on the dual formulation. 

This formulation has an advantage in that the solution space for 
zero penalty is greatly expanded even if the error range is small. This 
has improtant consequences for the implementation and achievement of 
acceptable solutions for problems with many constraints and variables. 
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VI. QUADRATIC PROGRAMMING FORMULATIONS 



A. GENERAL 

The quadratic programming model described in Grinold [Ref. 1] pro- 
vides an alternate method for optimization of allocations. A general 
discussion of the solution techniques for quadratic programs is found in 
Simmons [Ref. 9]. The following programs are expressed with only equality 
constraints and will yield solutions that are found by analytically 
solving systems of equations. This provides a useful and consistent 
means of determining allocations. However, the use of QP as presented 
here must be weighed against two disadvantages. 

• The penalty weights must be symmetric for both deficit and excess 
types of errors. 

• The decision variables are unbounded and, thus, the QP may yield 
impractical solutions such as negative accessions. 

Since much of the data input for the model is subjective or highly 
variable, and the purpose of the model is to explore questions concerning 
policy planning, the magnitudes and relative relationships of the alloca- 
tion are more important than an exact choice of allocation. Thus, the 
use of symmetric penalties may not present significant problems. If the 
parameters and data inputs to the model are fairly consistent with real 
world conditions, then the use of unbounded variables is unlikely to 
produce inconsistent solutions which would adversely affect the result of 
the analysis. The appearance of negative values in the solution is 
likely to indicate misspecif ication of parameters or the attempt to 
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achieve unrealistic goals. For the following sections, it will be 
assumed that t|) = 0 and 6 = ii> . 

The QP is similar to a least squares approach in that the purpose is 
to minimize the sum of squared weighted errors. Appendix B contains a 
brief summary of the resulting quadratic forms used in this section. 
Reindexing the variables makes it possible to conveniently express the QP 
as a system of equations using matrix arithmetic notation. For mathemat- 
ical convenience, the following change of variables may be necessary: 

m = J (i - 1) + j (8) 

where m=l,2,...M and M=IJ; 

n = K (i - 1) + k (9) 

where n=l,2,...,N and N=IK. 

In this way, matrices of size (I,K) may be reindexed to become vectors 
of size (N) and matrices of size (I,J) become vectors of size M. Three 
dimensional arrays may also be reduced to two dimensions. For example, a 
three dimensional array of size (I,J,K) indexed over (i,j,k) may be 
reindexed as a two dimensional array of size (M,K) using equation (8). 

In using matrix notation, it will be necessary on occasion to reform a 
vector into a diagonal matrix. This means that the k^" element of the 
vector will occupy the (k,k)^^ position of the matrix and all off 
diagonal elements will be zero. 

B. ACCESSION AND INVENTORIES 

The first formulation to be considered is QPl(Y). The advantage of 
this formulation is that it requires less data input since the job 
sharing fractions (F) and requirements by rank and manoower class (P) are 
not needed. Also, once the optimal accession and inventories are 
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determined, values for F and P may be computed based on the optimization 



results. The objective function is: 



min z z u^.j( z w-j^ - b-j ) 



w.|^ and then reindex (i,j) to (m) and (i,j,k) to 
(m,k) as given above in equation (8). This results in: 



Now let h^-ji^ 



^ki j 



2 

min z u ( z h, y. - b„) 
m , mk ■'k m' 
m k 

which, in matrix notation, equals: 

min (HY - B)' U'U (HY - B) 

where U is now an (M,M) diagonal matrix of the penalty weights. Taking 
the derivative with respect to Y and setting the result equal to zero 
yields the following solution: 

Y = (H'U'UH)"^ H'U'UB . 



Once an optimal value for Y has been determined, inventories may be 
calculated using equations (7) and (6) from Section II. 

The formulation QP3(Y) minimizes the objective: 
min £ z u2.( - b,j)2 * I z 

subject to the sharing constraint 



”10 ' I 9kij "ik A • 

Substituting the constraint into the objective, letting h.j|^ = gj^^j w^-|^, 
and then reindexing from (i,j,k) to (m,k) as given above in equation (8), 
yields in matrix notation: 

min (HY -B)' U'U (HY - B) + Y'QY - 2Y'R + c 
where Q is a (K,K) diagonal matrix whose diagonal elements are given by 

V = ' 4 4 ; (10) 
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2 

R is a (K) vector with elements f ; 

2 2 

and c is a scalar constant, namelv c = t 2 v.. p.. . 

i k 

Taking the derivative and setting it equal to zero results in 
H'U'UHY - H'U'UB + QY - R = 0 . 

Finally, the solution is: 

Y = [H'U'UH + Q]‘^ [H’U'UB + R] . 

This model can be found in Grinold [Ref. 1] as UQ-P. An advantage to 
the use of this model over other QP3 formulations is that it uses the 
people sharing rule G which may be more available as a data input than 
the job sharing rule F. 

An alternative program which uses the same objective function is 
QP3(X,Y). In this case, minimize the objective subject to the sharing 
constraint: 



“ik 'kij ’'ij = ° 

Both X and Y are now the decision variables and the constraint cannot 
be substituted into the objective function directly. In order to express 
the constraint in matrix notation, it is necessary to define two new 
matrices, F* and W*. First, F* is defined by the elements: 



* 

^i 'ki j 



^ij ^ 



0 otherwise 
where i, i' k=l,...,K, and j=l,...,J. 

Now reindex from (i',k,i,j) to (n,m) as given by equations (8) and (9) 
above to form the properly shaped F*. This results in a block matrix in 
which each of the blocks is a diagonal matrix of the j^*^ column of F 
for a particular (i,k). W* is formed in a similar manner. Let: 
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r 



for k = k 



0 otherwise 
where k, k'=l,...,K and 

Next reindex W* from (i,k',k) to (n,k) according to (9) above. 

The result is a diagonal block matrix in which each of the blocks on 
the diagonal is a column of W. Finally, the constraint can be expressed 
in matrix notation as: 

F*X - W*Y = 0 . 

The objective function in matrix notation is: 

min (X - B)' U'U (X - B) + Y'QY - 2Y'R + c 
where U, Q, R, and c are the same as defined for QP3(Y). Now form the 
Lagrangian: 

(X - B)' U'U (X - B) + Y'QY - 2Y'R + c - X (F*X - W*Y) 



where x is the vector of Lagrange multipliers. Taking the derivative 
with respect to X, Y, and x yields the following system of equations: 



2U'UX - 2U'UB - F*' 


= 0, 








2QY - 2R + 'W* 


= 0, 








F*X - W*Y 


= 0. 








This can be written in partitioned matrix form 


as: 
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where 0 indicates a zero matrix of appropriate size. The optimal Y, X, 



and \ is given by; 
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The program QP3(X,Y) is found in Grinold [Ref. 1] with the name UQ-B. As 
can be seen from the result, QP3(X,Y) has a more complicated solution 
than QP3(Y) and has an additional disadvantage in that the results may 
slightly violate the consistency rule. Thus, a determination must be 
made whether to use both X and Y from the results or to construct the 
inventories solely on the basis of Y. It has been presented here as an 
illustration of the use of Lagrange multipliers and as a means of checking 
the results from QP3(Y). 

Another means of checking QP3(Y) which insures that the consistency 
rule is met would be to solve QP3(X), use equation (5) to determine z, 
and then use the least squares approach presented in Section II to find 
values for Y. Those values of Y may then be used to reconstruct the 
inventories and insure consistency. 

C. EXPECTED STAGE LENGTHS 

The next formulation that will be discussed is QP2(W). The objective 
function is: 
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' 4 "ik - Pik)^ 

w 1 k 

Oik “ 4 ^k 0"0 '-ik ' "ik Pik • 

Then the result is: 

; (Oik “ik - '"ik>^- 

w 1 k 

Now re indexing from (i,k) to (n) and changing to matrix notation yields: 
min (QW - R)' (QW - R) 

where Q is reformed to be an (n,n) diagonal matrix. The solution to 
this program is: 

W = [Q'Q]"^ Q'R 

0. ADDING ADDITIONAL EQUALITY CONSTRAINTS 

The technique of Lagrange multipliers presented with the QP3(X,Y) 
model may also be used with additional equality constraints. For example, 
suppose the total number of accessions to be brought into the system is 
known but the exact number in each manpower class is not known. This 
information may be added as an additional equality constraint: 

2 Y;, = r 

where r is now the total number of recruits to be brought into the 
system each time period. In matrix notation, this constraint becomes: 
c'Y - r = 0 

where e is a vector of k ones. Using the QPl(Y) model as an example, 
form the Lagrangian with this additional constraint to get 
min (HY - B)' U'U (HY - B) - X (c'Y - r) . 
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Take the derivatives, set equal to zero, and solve as before to get 



the solution: 
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1 

L 


1 

J 



The advantage of adding an additional constraint is that it can be 
used to force the quadratic model to search for solutions that will be 
more related to "real world" problems. This helps to mitigate the 
disadvantage regarding unbounded decision variables. In fact, a range of 
possible values could be specified and then the optimal solution within 
the range could be found using a simple numerical search procedure. 

Since the solution to each problem during a search iteration is explicit 
and the resolution needed is to the nearest whole integer, such a proce- 
dure is easy to implement and should arrive at solutions fairly quickly. 
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VII. DEMONSTRATION OF THE MODEL 



A. IMPLEMENTATION 

The model was implemented in its quadratic form on an IBM 3033 
computer at the Naval Postgraduate School using the APL programming 
language as shown in Appendix C. For the linear and elastic programs, 
preprocessors were written in APL to arrange the data in standard MPS 
programming format for implementation on the XS-system optimization 
package. MPS is an international standard and such format is easily 
transferrable to other types of computers. A sample output is included 
in the program listing section. In addition, listings of some of the APL 
program functions used to produce the results in this section are also 
included. The ease with which data arrays can be manipulated in APL 
resulted in concise and flexible program functions which performed the 
optimizations and the input-output formatting of the data. Also written 
were two "user friendly" functions for data input and data display, and a 
master driver function that allows the user to select from a menu of 
options. Appendix C contains a brief summary of these programs. 

B. SYSTEM DESCRIPTION 

1 . Attributes 

This demonstration is based on an example from the U.S. Navy 
Officer Corps. The attributes chosen for manpower classes, career stages, 
and job classes are the same as those found in Grinold [Ref. 1] and are 
shown in Figure 13. 
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MANPOWER CLASSES 



Type 

1 

2 

3 

4 

5 



Officer Type 


Description 




IlOx 


General Unrestricted 


Line 


lllx/116x 


Surface Warfare 




112x/117x 


Submarine 


Warfare 




131x/139x 


Pilots 






132x/137x 


Naval Flight Officer 


(NFO) 




CAREER STAGE 






Stage 


Rank 


Years of 


Servi 


1 


ENS 


0 - 


2 


2 


LTJG 


2 - 


4 


3 


LT 


4 - 


9 


4 


LCDR 


9 - 


14 


5 


COR 


14 - 


19 


6 


CAPT 


19 - 


26 




JOB CLASSES 







Type Billet Designator 



Description 



1 1000 

2 1050 

3 1110/1160 

4 1120/1170 

5 1310/1390 

6 1320/1370 

7 1300 



General, Nonwarfare 
General , Warfare 
Surface Warfare 
Subsurface Warfare 
Pilots 

Naval Flight Officer 
General Aviation 



Fig. 13. Attributes for Example. 



2. Expected Time in Each Rank 

To compute W, it is necessary to have stage lengths and survivor 
fractions. Navy usage of survivor fractions is usually found in the form 
of continuation or retention rates. Current data was obtained from 
Deputy Chief of Naval Operations (MPT), OP-01, with continuation rates 
for Navy Unrestricted Line. The continuation rates were used directly to 
compute expected length of time spent in each rank. This data is shown in 
Appendix D. The APL function WAITS uses continuation rates and stage 
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lengths as inputs and produces the matrix W. the program GWAITS will 
also accomplish this but uses survivor fractions instead of retention 
rates. The resulting steady state matrix W is shown in Table 8. 



TABLE 8 

EXPECTED STAGE LENGTHS 





GURL 


SURF 


SUB 


PILOT 


NFO 


ENS 


1.925 


1.942 


1.891 


2.000 


1.992 


LTJG 


1.674 


1.700 


1.621 


1.975 


1.887 


LT 


2.063 


2.073 


2.040 


2.827 


3.241 


LCDR 


1.106 


1.138 


1.040 


1.028 


2.255 


CDR 


0.856 


0.896 


0.777 


0.773 


1.869 


CAPT 


0.596 


0.589 


0.607 


0.445 


1.276 



3. Billet Requirements and Sharing Rules 

Information concerning the requirements by rank and job type was 
also obtained from OP-01 and the array B was constructed. (See Table 
9). The next parameter to be specified is the billet sharing array. This 
information was taken from Grinold [Ref. 1] with some slight modifications 
and is displayed in Table 10. Once B and F are available, it is possible 



TABLE 9 



BILLET REQUIREMENTS BY RANK AND JOB TYPE 





1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


330 


0 


2258 


660 


1105 


593 


0 


LTJG 


768 


0 


1706 


700 


1970 


1185 


2 


LT 


1780 


378 


2145 


884 


3780 


1408 


629 


LCDR 


1456 


571 


1510 


877 


1818 


586 


716 


CDR 


1001 


468 


990 


518 


710 


54 


909 


CAPT 


571 


462 


354 


152 


0 


0 


422 



54 







I 



( 



i 



TABLE 10 



BILLET SHARING ARRAY 



GURL 


1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.73 


0 


0 


0 


0 


0 


0 


LTJG 


0.70 


0 


0 


0 


0 


0 


0 


LT 


0.70 


0 


0 


0 


0 


0 


0 


LCDR 


0.70 


0 


0 


0 


0 


0 


0 


COR 


0.70 


0 


0 


0 


0 


0 


0 


CART 


0.70 


0 


0 


0 


0 


0 


0 


SURF 


1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.18 


1.00 


1 


0 


0 


0 


0 


LTJG 


0.15 


0.50 


1 


0 


0 


0 


0 


LT 


0.15 


0.50 


1 


0 


0 


0 


0 


LCDR 


0.15 


0.50 


1 


0 


0 


0 


0 


CDR 


0.15 


0.50 


1 


0 


0 


0 


0 


CART 


0.15 


0.50 


1 


0 


0 


0 


0 


SUB 


1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.09 


0.00 


0 


1 


0 


0 


0 


LTJG 


0.06 


0.10 


0 


1 


0 


0 


0 


LT 


0.06 


0.20 


0 


1 


0 


0 


0 


LCDR 


0.06 


0.20 


0 


1 


0 


0 


0 


CDR 


0.06 


0.20 


0 


1 


0 


0 


0 


CART 


0.06 


0.20 


0 


1 


0 


0 


0 


RILOT 


1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.00 


0.00 


0 


0 


1 


0 


0.00 


LTJG 


0.50 


0.10 


0 


0 


1 


0 


0.00 


LT 


0.50 


0.18 


0 


0 


1 


0 


0.57 


LCDR 


0.50 


0.18 


0 


0 


1 


0 


0.57 


CDR 


0.50 


0.18 


0 


0 


1 


0 


0.57 


CART 


0.50 


0.18 


0 


0 


1 


0 


0.57 


NFO 


1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.00 


0.00 


0 


0 


0 


1 


1.00 


LTJG 


0.40 


0.10 


0 


0 


0 


1 


1.00 


LT 


0.40 


0.12 


0 


0 


0 


1 


0.43 


LCDR 


0.40 


0.12 


0 


0 


0 


1 


0.43 


CDR 


0.40 


0.12 


0 


0 


0 


1 


0.43 


CART 


0.40 


0.12 


0 


0 


0 


1 


0.43 



55 



to calculate the corresponding requirements in terms of people types and 
ranks, and the corresponding people-sharing array. The results for the 
requirements by people types is shown in Table 11. An example of the 
people-sharing array is shown in Table 12 for the manpower class of 
pilots. 



TABLE 11 



BILLET REQUIREMENTS BY RANK AND MANPOWER CLASS 





GURL 


SURF 


SUB 


PILOT 


NFO 


ENS 


241 


2317 


690 


1105 


593 


LTJG 


538 


1821 


746 


2010 


1217 


LT 


1246 


2601 


1066 


4296 


1795 


LCDR 


1019 


2014 


1079 


2402 


1021 


CDR 


701 


1374 


672 


1362 


541 


CAPT 


400 


671 


279 


352 


260 



TABLE 12 

PEOPLE-SHARING ARRAY FOR PILOTS 





1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


0.015 


0.000 


0 


0 


0.985 


0 


0.000 


LTJG 


0.019 


0.000 


0 


0 


0.980 


0 


0.001 


LT 


0.021 


0.016 


0 


0 


0.880 


0 


0.083 


LCDR 


0.030 


0.043 


0 


0 


0.757 


0 


0.170 


CDR 


0.037 


0.062 


0 


0 


0.521 


0 


0.380 


CAPT 


0.081 


0.236 


0 


0 


0.000 


0 


0.683 



4. Permitted Unit Errors 

The permitted unit errors were determined based on subjective 
judgment. First the matrix for permitted unit errors by rank and job 
type was estimated. This is shown in Table 13. The entries marked as 1 
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TABLE 13 



PERMITTED UNIT ERRORS BY JOB TYPE AND RANK 





1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


400 


1 


500 


500 


500 


400 


1 


LTJG 


400 


1 


500 


500 


500 


400 


500 


LT 


10 


8 


8 


10 


8 


10 


10 


LCDR 


10 


8 


10 


15 


8 


8 


10 


CDR 


10 


8 


5 


10 


8 


10 


10 


CAPT 


10 


10 


10 


10 


1 


1 


10 


correspond to 


ranks and 


job 


types for 


which 


there are 


no billets. 


Giving a very 


low percentage 


will help 


ensure the optimization 


does not 


seek to place 


people in 


these jobs. The corresponding 


errors 


by people 



type and rank were generated by the use of the consistency rule. The 
unit errors do not necessarily have to meet this requirement, but it was 
felt that the consistency rule would give an approximation as reasonable 
as further use of subjective judgment. The results were checked to insure 
that the numbers were reasonable, and no changes were made. The matrix 
for permitted unit errors by rank and people type is shown in Table 14. 
These last two matrices represent the permitted error for overfilling the 
billets. The matrices for underfilling the billets were done using a 
similar method. 

C. ACCESSIONS USING QUADRATIC PROGRAMMING 

The first set of programs to be considered are those labelled QPl(Y), 
QP2(Y), and QP3(X,Y). The results of the optimization with respect to 
the accessions are shown in Table 15. The column marked "Obj . Val." is 
the value of the penalty function that measures both types of errors. 
Notice that QPl(Y) which only minimized with respect to errors by job 
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TABLE 14 



PERMITTED UNIT ERRORS BY RANK AND PEOPLE TYPE 





GURL 


SURF 


SUB 


PILOT 


NFO 


ENS 


280.0 


560.5 


524.2 


520.7 


416.5 


LTJG 


280.0 


560.5 


524.2 


805.2 


631.1 


LT 


7.0 


13.5 


12.2 


15.6 


15.7 


LCDR 


7.0 


15.5 


17.2 


15.6 


13.7 


CDR 


7.0 


10.5 


12.2 


15.6 


15.7 


CAPT 


7.0 


16.5 


12.6 


9.0 


6.9 



TABLE 15 

ACCESSION USING QUADRATIC PROGRAMMING 
Obj. 





Val . 


GURL 


SURF 


SUB 


PILOT 


NFO 


Total 


QPl(Y) 


455.2 1 


i 704 


1404 


554 


1470 


332 


1 4464 


QP3(Y) 


390.6 1 


i 712 


1416 


563 


1246 


272 


1 4209 


QP3(Y,X) 


400.8 1 


i 704 


1380 


556 


1388 


270 


1 4298 



type and rank still did well when compared to the other two. In fact, as 
can be seen from the actual accession values, all three programs had 
remarkably close results. Since QP3(Y) did the best in terms of the value 
of the objective function, the results from this optimization will be 
presented in the remainder of this section. 

While close, the results from all three programs have rather high 
penalty values. Recall that the minimum value that can be achieved 
is zero, and the actual magnitude of the penalty values shows that our 
actual allocation is far off the mark. Another way to gain an apprecia- 
tion of how close the allocation is to the targets is by examining the 
percent error compared to the permitted errors. Listed below in Table 16 
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TABLE 16 



PERCENT ERROR BY RANK AND JOB TYPE 





1000 


1050 


1110 


1120 


1310 


1320 


1300 


ENS 


351 


0 


19 


54 


125 


-9 


0 


LTJ6 


90 


0 


32 


22 


22 


-58 


-12 


LT 


12 


-1 


13 


8 


-18 


-51 


-32 


LCOR 


-26 


-32 


-20 


-46 


-47 


-40 


-44 


CDR 


-14 


-17 


-8 


-35 


-29 


-6 


-19 


CAPT 


14 


31 


24 


23 


0 


0 


47 


are the percent errors 


from the QP3(Y) 


program. 


This 


table highlights 



problem that is probably due to the high rate of attrition of the middle 
level officers. Notice that those ranks are underfilled (as indicated by 
the minus sign) while the junior and senior ranks are experiencing an 
overfill. The most serious problem exists in the rank of Lieutenant 
Commander and in the job classes Submarine Warfare Officer and Pilot. 

The model could be made to actually fill those shortage billets by moving 
the permitted unit error in those ranks and job classes closer to zero. 
However, since the current parameters are already closer to zero than any 
of the others, such changes would only result in worse overfills at the 
high and low end. A better approach would be to revise the sharing rules 
in an attempt to get a better allocation. 
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VIII. SUMMARY 



A. GENERAL 

The purpose of this thesis was to investigate and enhance the model 
proposed by Grinold [Ref. 1] with special attention to the optimization 
techniques and considerations involved with implementation. The model 
looks at requirements and inventories of a manpower system by classifying 
them according to people type, rank, and job type. By the use of sharing 
fractions, which represent management policy, it is possible to express 
billets which are usually defined by job type as requirements for certain 
types of people. The same sharing rules also make it possible to express 
inventories of people by the types of jobs filled rather than by the 
types of people filling those jobs. Further, the assumption of steady 
state makes it possible to calculate all inventories from accessions. 

Thus, billets, inventories, and accessions are linked by a unified and 
rather elegant model structure. 

If relative trade-offs between different billets in the form of 
permitted shortfalls or overfills can be made, then it is possible to 
measure the error between the desired billet structure and the steady 
state inventories. Various optimization techniques may then be used to 
examine the long term result of policy decisions regarding the allocation 
of people to jobs and to test alternative policies. 

This thesis has included several areas of research not presently 
found in Grinold [Ref. 1]. Several considerations were discussed concern- 
ing the implementation of the model. These include: level of detail in 
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attribute definitions, problems encountered in implementation, acquisition 
of model parameters, intuitive explanation of weighting factors, and 
advantages of various optimization schemes. New programming formulations 
were presented that used accessions, inventories, sharing rules, or 
expected stage lengths as the decision variables. Elastic programming 
was introduced as a generalization of the goal programming technique. 
Explicit solutions were derived for the quadratic programming model and 
the idea of additional equality constraints was proposed. Other enhance- 
ment features included the least squares approach to determine the 
requirements by people type and the job sharing rules given the billet 
structure and the people sharing rules, and using least squares to 
determine accessions from a given inventory. 

In conclusion, this model allows the policy planner to examine a wide 
range of policy options using only a few simplifying assumptions and a 
modest amount of input data. Many areas of policy planning may be 
illuminated by simply placing a problem in the framework of the model 
since this requires that the analyst consider trade-offs between various 
people and job types and formalize policy in setting the values for the 
sharing fractions. This ability to examine the aggregate effect of human 
behavior on an organization's billet structure makes this model a useful 
tool for the policy planner. 

B. AREAS FOR FURTHER STUDY 

There are several areas related to this model where further research 
should be done. Suggested areas include the following; 

• Extend the model structure to handle dynamic conditions. If 
the current inventory is known and the continuation rates of 
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personnel can be estimated, then a year-by-year approach could 
be used to examine various short-range policy problems. This 
would eliminate the limitation of steady state assumptions in the 
current model . 

• The optimizations rely heavily on reasonable estimates of the 
permitted percent error in failing to fill all billets. More 
research is needed on how to gather, scale, and interpret data 
that could be used by the model with regard to these permitted 
errors . 

• More improvement is needed in adding cost considerations to the 
model. The ability to relate any analysis to "the bottom line" 
should enhance the results of such an analysis. The idea of a 
multi-attribute objective function that considers both errors and 
costs was presented in Section III. Another possibility would be 
to develop a means to econometrical ly express the trade-off 
parameters. 

• The usefulness of the model would be extended if the results could 
be expressed in such a manner that the impact of policy decisions 
could be quickly evaluated. One idea is to use the model to 
create ratios or indices that pinpoint the criticality of various 
allocations. For example, if the "optimal" billet structure based 
on projected inventories were calculated, then the ratio of this 
theoretical structure and the actual structure should provide a 
measure of "how critical" certain billets will be in the future. 
Ratios close to one are desirable, and those far away hignlight 
the problem areas. 
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• Since the model is based on steady state assumptions, it is 
desirable to be able to relate the results of the model to these 
assumptions. Assuming certain functions that describe the way in 
which people leave the organization could be used to establish 
such a relationship. The impact of changes in continuation rates 
or survivor fractions on the billet structure needs to be directly 
related to the sharing policies that will be required to reduce 
the allocation errors, and conversely, changes in the billet 
structure need to be analyzed in light of the behavior of people 
in the organization. 

• Further model enhancement is needed in order to introduce the 
effects of stage substitution to the model. In military organiza- 
tions, such stage substitutions exist in the form of "selected but 
not yet promoted." The model has assumed that such stage substitu- 
tion (i.e., personnel of one rank filling a job of a higher rank) 
is negligible or has a net effect of zero. The ability to model 
this would add more credibility to the model and subsequent 
analysis. 
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APPENDIX A 





VARIABLE LIST FOR THE MODEL AND PROGRAMS 


Name Size 


Description 


k 


The index of the manpower classes. k=l,2,...,K. 


i 


The index of ranks or career stages. i=l,2,...,I. 


j 


The index of the job classes. j=l,2,...,J. 


B I,J 


The number of people in rank i required to fill class j jobs. 


P I,K 


The number of required people in rank i and manpower class k. 


T K,I,J 


The number of (i,k) people required in type (i,j) jobs. 


1 

X 


The inventory of people in rank i who are filling job j. 


Z I,K 


The inventory of people in rank i and manpower class k. 


A K,I,J 


The inventory of (i,k) people who fill type j jobs. 


F K,I,J 


The fraction of (i,j) jobs that should be performed by 
manpower class k. 


G K,I,J 


The fraction of type (i,k) manpower who should fill job 
type j. 


W I,K 


The expected length of time a person in manpower class k 
will spend in rank i. 


S I,K 


The fraction of type k people who will remain in the organ- 
ization past rank i. 




The permitted unit error for underfilling job class j 
and rank i . 


9 I,J 


The permitted unit error for overfilling job class j and 
rank i . 


« I,K 


The permitted unit error for underfilling rank i and man- 
power class k. 


4> I,K 


The permitted unit error for overfilling the rank i and 
manpower class k. 




The lower bound of permitted error for job class j and 
rank i . 
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B 


1,0 


The upper bound of permitted 
rank i . 


error for job class j and 




£ 


I,K 


The lower bound of permitted 
rank i . 


error for manpower class k 


and 


■p 


I,K 


The upper bound of permitted 
rank i . 


error for manpower class k 


and 


C 


K , 0 


An indicator that is 1 if manpower class k can fill job 
0 otherwise. 


j; 


e 


1,0 


^ U 

The error in meeting (i,j) 


requirement. 




d 


I,K 


The error in meeting (i,k) " 


requirement . 




u 


1,0 


The weighting factor for an ( 


i,j) error. 




V 


I,K 


The weighting factor 'for an ( 


i,k) error. 
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APPENDIX B 



QUADRATIC FORMS FOR THE MODEL 



1 . The General Form 

The quadratic objective functions in this thesis have the same general 
quadratic form which can be expressed in summation notation or matrix 
notation, if properly reindexed. Consider the quadratic form with 
variables X(I,1), A(I,1), P(I,1) and constant r: 



where Q is the diagonal matrix: 




2 . The Special Form 

If the above general quadratic form has the following characteristics, 
then a more familiar form can be written for it. 




(al) 



In matrix form this becomes; 



X'QX - 2X'P + r 



(a2) 




and then (al) becomes; 





(a3) 
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In matrix notation this is: 

X’D'DX - 2X’D'B + B'B 
(DX - B)’ (DX - B) 
where 0 is the diagonal matrix: 



D 



1 






‘I -1 



This special form has several cases which are of interest. 
Case 1 

Let a^ = 1. Then (a3) becomes: 



z (x . - b . ) ' 



(a4) 



(a5) 



or in matrix form: 

(X - B)’ (X - B) 

Case 2 

Let b^. = a^.c^ . Then (a3) becomes; 



z a^ (x. - c.)^ 
i 



X'D'DX - 2X'D'DC - C'D'DC 
(DX - DB)’ (DX - DB) 

(X - O' D'D (X - C) . 



Case 3 



(a6) 



(a7) 



(a8) 



Let b, = and x, = . 

Since FZ = X, this form is like (a7) and has similar results: 
£ a .2 [(£ - c ,]2 

or in matrix notation; 

(FZ - O' D'D (FZ - C) . 



(a9) 



(alO) 
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APPENDIX C 



APL FUNCTIONS 



Function Uses Computes 

QUAD All variables 



QP3XY B, F, P, W, Y, X, Z, A 

9, (fr 



MAKEPG P, G B, F, T 

MAKEBF B, F P, G, T 

MAKEFP B, G P, F, T 

QP3Y G, B, F, P, Y, X, Z, A 

W, 9 , <1) 

INVENTORY G, Y, W Z, X, A 



Syntax and Description 

QUAD. 

Master Driver Program. This 
function calls all others through 
the use of menu options. 

QP3XY. 

Calculates "optimal" accessions 
Y, personnel inventory Z, and 
allocation A, given the data, 
objective, and the billet share 
rule F. 

MAKEPG. 

Reconciles the global variables 
B, F, and T with the global 
variables P and G. 

MAKEBF. 

Reconciles the global variables 
P, G, and T with the global 
variables B and F. 

MAKEFP. 

Reconciles the global variables 
P, F, and T with the global 
variables B and G using least 
squares approach. 

QP3Y. 

Calculates "optimal" accessions 
Y, personnel inventory Z, billet 
staffing X, and allocation A, 
given the data, and use of the 
people share rule. 

INVENTORY. 

Calculates inventory Z, X and 
allocation A, where accessions, 
people share rule, and expected 
stage lengths are given. 



100(D-Z)/P 



PEOPRCNT P, Z 
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PEOPRCNT. 

Calculates the percentage error 
in inventory. Compare with 4>. 




( 



Function Uses 



Computes 



Syntax and Description 



BILPRCNT 

GWAITS 

WAITS 

ROUND 

LABIJ 

LABIK 

LABKIJ 

DISPIN 

DISPOUT 



B, X 100(B - x)/B 



ALPHA, S W 



RET, S W 



ILAB, JLAB, 
OP, CS 



ILAB, KLAB, 
DP, CS 



ILAB, KLAB, 
JLAB, DP, CS 



G, F, B, P, 
ILAB, KLAB, 
JLAB, OP, CS 



X, Y, Z, 8, 
ILABS, KLAB, 
JLAB, DP, CS 



BILPRCNT. 

Calculates the percentage error 
in meeting billet requirements, 
compare with 0. 

GWAITS. 

Given the survivor fractions 
ALPHA and stage definitions S, 
calculates W the expected waiting 
time in each stage. 

WAITS. 

Given the retention rates RET 
and stage definitions S, cal- 
culates W, the expected waiting 
time in each stage. 

ROUND XXX. 

Takes any array and rounds 
elements to integers. 

LABIJ XXX. 

Takes any array of dimension 
(I,J) and labels rows and col- 
umns according to ILAB an JLAB. 
ILAB must have I rows and JLAB 
must have J rows. OP specifies 
number of decimal places desired. 
CS specifies minimum column 
spacing. 

LABIK XXX. 

Handles any array of dimension 
(I,K) in the same manner as LABIJ. 

LABKIJ XXX. 

Handles any array of dimension 
(K, I, J) in the same manner as 
LABIJ. 

DISPIN. 

User selects from a menu to 
produce the input variable with 
appropriate labels. 

OISPOUT 

User selects from a menu to 
produce the output variable with 
appropriate labels. Also displays 
error arrays from 3ILPR0NT and 
PEOPRCNT. 
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APPENDIX D 



NAVY URL CONTINUATION RATES ^ 



Years 


GURL 


SURF 


SUB 


PILOT 


NFO 


1 


0.963 


0.971 


0.947 


1.000 


1.000 


2 


0.960 


0.971 


0.938 


1.000 


0.985 


3 


0.930 


0.928 


0.933 


0.990 


0.970 


4 


0.820 


0.809 


0.841 


0.980 


0.919 


5 


0.698 


0.698 


0.697 


0.882 


0.858 


6 


0.816 


0.819 


0.810 


0.692 


0.882 


7 


0.870 


0.860 


0.891 


0.706 


0.859 


8 


0.912 


0.914 


0.908 


0.808 


0.959 


9 


0.947 


0.970 


0.902 


0.818 


0.972 


10 


0.843 


0.833 


0.862 


0.828 


0.934 


11 


0.868 


0.884 


0.835 


0.883 


0.918 


12 


0.922 


0.922 


0.922 


0.933 


0.927 


13 


0.914 


0.926 


0.890 


0.969 


0.971 


14 


0.981 


0.987 


0.970 


0.964 


0.986 


15 


0.974 


0.971 


0.979 


0.954 


0.974 


16 


0.989 


0.990 


0.987 


0.948 


0.974 


17 


0.978 


0.974 


0.987 


0.957 


0.970 


18 


0.956 


0.964 


0.941 


0.950 


0.936 


19 


0.917 


0.901 


0.948 


0.838 


0.954 


20 


0.762 


0.754 


0.778 


0.765 


0.789 


21 


0.831 


0.821 


0.851 


0.774 


0.692 


22 


0.908 


0.883 


0.957 


0.901 


0.914 


23 


0.875 


0.853 


0.920 


0.844 


0.860 


24 


0.911 


0.894 


0.946 


0.852 


1.000 


25 


0.848 


0.847 


0.849 


0.838 


1.000 


26 


0.701 


0.710 


0.682 


0.778 


1.000 



^Source: DCNO, OP-01. 
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APPENDIX E 



COMPUTER PROGRAMS 



VBDISPLAYlUl^ 

7 dDISPLAYiS 

'EUTER SYMBOL OF VARIABLE TO BE DISPLAYED' 

' TO QUIT TYPE 0* 

f f 

’ 3 (.BILLET REQUIREMENTS) ' 

' F (JOB SHARING FRACTIONS) ' 

' T (TARGET ALLOCATIONS)' 

' X (CALCULATED BILLETS FILLED)' 

' A (CALCULATED ALLOCATIONS) ' 

' E (PER CENT ERROR IN MEETING BILLET REQUIREMENTS) ' 

L0:*( (S^a) = ' BFTXAE' ) / L-i , L2 , L2 , L^ , L2 , L2 
-*■0 

Llzt'LABIJ ' ,S 
■*L0 

L2:t' LABKIJ ' ,5 
•*L 0 

L3i I'LABIJ ROUND BILPRCNT 
■*L 0 
7 

\;BILPRCNTiWi^ 

7 BILPRCNT 
LI] LABIJ ROUND 1 OQx (£-3 ) tS+3= 0 
7 



7SJiV$PCQ]7 
7 Z*-3INQP RTOT 

l1 ] H*-( ( 1 ) + . . x3 1 ) ,-0.5x£pl 

L2] *T^P,[l](£pl ) ,0 

l 3] 7i:^( ( (S?31) + .x^+.x ,3) ,[i ] RT0T)'SH 

L 4 ] ^ 

L5] INVENTORY 
[ 6 ] Z^PF\ 

7 



VBINSEARCHiWiN 
7 X^BINSEARCH P 
Ll] £C(?P;-(gP5>|P[l ]-?[23)/0 
L2] A'-^( + /P)t2 

13] ^(Q<3INQP a)/£1 

14] P^2'\^X 

15] '*L00P 

16] u1:PL1] '*'X 
L7] ^LOOP 

7 



[ 1 ] 

C2] 

L3] 

L4] 

L5] 

[ 6 ] 

C7j 

L8] 

L9] 

Lio] 

Cll] 

l12] 

LI 3] 

[ 14 ] 

[ 15 ] 

[16] 
l17] 
[18] 
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'^CHAi! C/cC u ] V 

V CHAUGE 

l1] 'y^HAT IS THE uA:^3 OF THE VARIABLE TO BE CHARGED? • 

L 2 ] 

L3] 0* ) /END 

l4] ^((p,A)= 1 3 5) /Cl ,(72 , (73 

L5] ■*ERN2 

16] C2ttIF *x/A'= * * PHI * ♦ * THEN ’ ’ ELSE '-*ERR2' 

17] C3:tIF *x/A'=* ' THETA ' ’ ' THEN '^L' ELSE ^^ERR2' 

LU] 3PFGW ' ) /ERR2 

[9] LiWHAT IS THE LOCATION OF THE ELEMENT YOU NISH TO CHANGE?^ 
LlO] i 0;-^(0 = x/(t»p,* ,X)^KIJ^Q) /ERRr 
ill] 'INPUT NEN ’ ,A, ’ ( ’ , (f AJJ’) , ’ ) ’ 

[ 12 ] iNDEX^'fXIJil^ll 

[13] L1i^{(.qKIJ)<I^I^\) /OUT 

[14] INDEX^IUDEX i' ,^KUiI'^ 

[15] -Cl 

[16] OUT:tX,'l' .INDEX.'I^U' 

[17] -*END 

[18] ERR^ INPUT ERROR. TRY AGAIN.' 

[19] -CO 

[20] ERR2'.' SYMBOL NOT UNDERSTOOD' 

[21] 'ENTER AGAIN OR TYPE 0 TO QUIT' 

[ 22 ] -1 

[23] END:' END OF CHANGE PROGRAM' 

7 



7C^//4C[G]7 
7 Z^CDIAG X\IJ\J 
L 1 ] IJ-*-qX 

[2] Z^DG A[;J’-1] 

[3] Cl :-( Jt/[2]<t/-./+l) /O 

l4] Z-Z,[1] DG XliJi 

[5] -Cl 

7 



7Cyi=l/[Ll]7 
7 CiPl/ 

[1] 'INPUT TOTAL NUMBER OF RECRUITS DESIRED' 

[2] RTOT^‘C\ 

[3] .¥-( ((7^1 ) + .x^+.x£l) ,-0.5xAp1 

[4] f/— , [ 1 ] ( Ap 1 ) , 0 

l 5] /C-( ( (i7tfl) + .x^+.x,5) ,[1 ] RTOT)m 

lS] l^X^YL 
7 
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ill 

C2] 

13] 

14] 



7 RR*- 'OG 

/V’-^P A 

* i/ ) P 0 
ifit'L 1 ; l-^-^ 

ixa *’ (3 ,“ii/"l)'l^ill 
7 



Rii 



Lll 

[2] 

C3l 

l41 

i5l 

Co] 

C7l 

C8l 

[9] 

Ciol 

-ill] 

[121 

Cl3l 

[14] 

[15] 

[16] 
[17] 
[13] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 

[25] 

[26] 
[27] 
[23] 



[ 1 ] 



7£?I5PI;V[lj]7 

e%Ir^ 5 ni BO L OF VARIABLE TO BE DISPLAYED* 
TO QUIT TYPE 0’ - 

TO SEE MEtiU TYPE M* 






f 

3 {BILLET REQUIREMENTS)* 

F {JOB SNARING ARRAY)* 

T {TARGET ALLOCATIONS)* 

G {PEOPLE SHARING ARRAY)* 

~P {PEOPLE REQUIREMENTS)* 

V- {EXPECTED FURTHER DURATION)* 



*13’ 

-*•(1 ^pS'*'t3) / ERRy 

•*{S= * BFTGPWM * ) /LI ,1.2 

•*END 

Ll:i*LABlJ ROUND ’ ,S 
*LQ 

L2:t*LABKIJ * ,S 
-LO 

L3ii*LABIK * ,S 



L5,L2,L6,L3,L4 



-LO 

LB’.i* LABKIJ ROUND ',5 
*LQ 

L6 :t* LABIK ROUND * ,5 
■^LO 

ERR^ : * INCORRECT ENTRY* 

END:* END OF DISPLAY PROGRAM* 



7 



7r^£;.v[Ll]7 

VIFCGIV ^ Z99-^/i99 THEN 599 

V Z^IF A 299-^,3991 1 +1^499 ; ] 

Z^A y 

7 



7LL53CG1V 
7 Z-^A ELSE 3;F 
Ll] F*-{p ,A)[ 9 tB 

c2] Z^(3+3) , [0 . 5] Ft/i 

7 
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: m 



1 



VG,r-AITSlQl7 

7 A I TS \A iSF iiJ ;J \AA\3B\DD\RH 

Cl] /i-^S 
L 2 ] 3'F^ ALFA 
il 3 ] i]'*- ^ fp5F 

[4] 

[5] AA-*-A « ,^in' 

L 6 ] B B-^-A » . S 1 + 1 // 
l 7 ] 33-*-B Bf^'^AA 

Cd] DD^^df ,I)pA-lA 
L9] AA'*~hA + DD^33 
LlO] DD->-{I ,1 ,JpO 

Lll] AA-*-DD*.xAA 

[12] RR*AA*.^SF 
L 1 3 ] ^-^RR 

7 



VIGVERTCRYiUl^ 

7 li/VERTORY 

[1] Z-^Wx(1,K)qY 

[2] /«^+/[l] ^«^£x«)(£,J.£)pZ 
7 



'JLABIJLOl'J 

7 Z^LABIJ XiN ilJ il iJ iL iH L;!'i \PP ;iJD iFS ;H iX,'J iXX iRS 

[1] flrA’JS PROGRAM LABELS Ail (I^J) ARRAY USIllG JLAB AtlD ILAB AS 

[2] nTHE LABELS FOR ROWS AtID COLUMi/S RESPECTIVELY . DP_ IS THE 
L3] rUUMBER of decimal PLACES TO BE DISPLAYED ADD C5 IS THE 

[4] rMIIUMUM SPACniG BETWEEU COLUMNS. 

[5] N*Q 

[ 6 ] I^(pa)[ 1 ] 

[7] c7-^(pX)[2] 

[8] ^(J«(pJi£B)[l])/5.7i?1 

[9] *(J2(p£££3)Cl])/Ei?i?2 

[10] - L^{p ILAB )i2'[ 

[ 11 ] HL-^i. o JLAB ){.2'^ 

[12] M^^U2,L)q' ' ) .[1 ] ,[2]( (J + 2) ,2)p’ |’ 

[13] £1 :-^(c7<//'^;7+l ) /i2 

[14] PP^D£^Ox+/i\XX^\Xiiin 

[15] Zv'i?^pT[//[( U/)[l ]] 

l 16] //5«^0^ + /( 0>;r[ ;^V] ) 

[17] FS^CS*HLd'JD^PP^NS-^l 

[18] /z/^(F5,?P)t(I,1)pA'L ;;/] 

[19] H^UFS-HDq’ • ) , JLAB lN',l 

[20] H^(2,pH)pH AqH)p' 

[21] M^M . [ 2 ] // , [ 1 ] XN 

[22] ^il 

l23] L2:Z*-M 
L 24 ] -*-0 

[25] ERR^ ERROR: NUMBER OF ROWS = NUMBER OF ROW LABELS' 

[26] -*-0 

l27] ERR2:' ERROR: NUMBER OF COLUMNS * NUMBER OF COLUMN LABELS' 

[23] -*-0 
V 
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L-J 



V 56 '- 4 ;-/LIJ] V 

V SSAM 

L 1 ] iJJTl MSU!/ D 0 ^ 
2 ] ' £HD OF SSAl-r 

V 



^ FSTAii LOl^ 

V F ^FSTAR \I 

Cl] 

L 2 ] /: 0 :^(J<I^I + 1 ) /'iil 

L 3 ] "*0j 0,p^C-^» • ]'*'^C t F t 1 

[ 4 ] £,1 (Jx^) . (I*£) 

V 



vPoCipifCi/rCGDv 
V PEOPliCUT 

[ 1 ] iS -4 5 /A: PCi/.VZ? 1 00 * (_Z-P) ^P+P =0 
7 



7 PP 1 lG]V 
7 P-*-PPl 

Ll] P^+/ + /(_i/x£-P )*2 

7 



7 PP 2 lU ]7 
7 P-^PF 2 

Ll] P^+/ + /(/x_ 2 -P )*2 
7 



7 PP 3 LU ]7 
7 ED*PF 2 

[ 1 ] dP^PP 1 +PF 2 
7 



7 i,'i>';iP[U ]7 
7 QUAD 

l1] t POQUAD l TYPE j ] 
C 2 ] I/,'/F// 2 ’PPr 

7 



7 /t£'i 7 iVPLu ]7 
7 nh^tiCUUD A A 

Cl] RR^{lAA) + {AA-lAA)iQ .5 

7 
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I 



Cl] 

l2] 

L3] 

L4] 

C5] 

L6] 

L7] 

L8] 

C9] 

LlO] 

[11] 

[ 12 ] 

[13] 

[14] 
[3 5] 
[16] 
[17] 
[13] 
[39] 
[ 20 ] 
[ 21 ] 
[ 22 ] 
[23] 
t24] 

[25] 

[26] 
[27] 
[23] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 

[35] 
l36] 
[37] 
[33] 

[39] 

[40] 

[ 41 ] 

[42] 

[43] 

[44] 

[45] 
L 46 ] 

[47] 

[48] 

[49] 
L 50] 
[51] 



VJ,VPOT[ Li] 7 
V Ti}PUTiKiKIJ;NiS;Z 

'Eii'IER THE HUMBER OF MAUPOWSR TYPES, HUMBER OF HARKS, 
'HUMBER OF BILLET TYPES AS A VECTOR.' 

f f 

'TYPE THE SYMBOL OF THE VARIABLE YOU VISH TO INPUT' 

' TYPE 0 TO QUIT' 

' TYPE M TO SEE MENU AGAIN' 

L:' 

' B_ BILLET REQUIREMENTS (JxJ)' 

'P PEOPLE REQUIREMENTS (I^K)' 

’£ BILLET SHARING ARRAY {K^I^J)' 

'G PEOPLE SHARING ARRAY iK^I^J)' 

'S SURVIVOR FRACTIONS U^K)' 

'R RETENTION RATES {N^K)' 

•e UNIT ERRORS IN TERMS OF BILLETS (.I^J)' 

’<}) UNIT ERRORS IN TERMS OF PEOPLE (.I^K)' 

f f 

lO : * ’ 

f f 



•n* 

-( Kp ,S^a) /ERR-i 
-(S*’ a » ) /LOl 
S^' THETA ' 

■*L 3 

LOl :-^(S«’4>* ) /LQ2 
S*- ' PHI ' 

-4 3 

L02:*iS='M BPFGSR ' ) / L , LU LZ , L2 , L2 , L^i , L^* 
■*END 

LI ‘.'INPUT EACH ROW OF ' ,S 
±S,'^MATIN KIJ 12 3]' 

-/:o 

:.2;;^-o 

Z^KIJpQ 

42 3 :-(a'JJ-[ 3 ]<;^-'^+1 ) /L22 
'FOR RANK ' ,fK 
'INPUT EACH ROW OF ' ,S 
ZIK‘, ',^^MATIN KIJ 12 3] 

-L21 

L22: tS ,'^Z' 

-*-L 0 

LZ:' INPUT EACH ROW OF ' ,S 
s.S,'^MATIN ^KIJi\ 2]* 

■*L 0 

L^’.' INPUT TOTAL NUMBER OF YEARS N' 

//-□ 

'FOR EACH YEAR, INPUT EACH ROW OF ' ,S 
iS,'^MATIH N,KIJi\'\' 

—4 0 

ERR^‘.' ERROR IN ENTRY OF VARIABLE' 

END'.' END OF INPUT PROGRAM' 

V 



AND' 
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I 



Vi'Pl 7 CJ ]7 

V 

l 1] ) + .*£^+.>‘.£)li](4?£1 ) + .x£y+.x^1 

V 



7iP3/CG]7 

V ^P3’/ 

Cl] Y^{R+{<SiHl)* .^LiU* ,5)a£+('S?^1 ) + .x££+. x^] 

V 



7^p3/rca]7 

7 QP3XY;MiRiH;XYC 

Cl ] ,K)pQ) ,-Q .5^<^ FSTAR 

C2] tf-^//.Cl]( (^..'<f)pO) .g.o.sx^ ^sr/ii? 

C3] .H:\i- FSTAR ) . WSTAR , (. R K) pO 

l4] XYL^i (.UjU* ,B) ,R,RpO)^H 
[5] Y_-*-M*{-(J)^XYL 

7 



7$P2t7CU]7 
7 QP2^ 

Cl] (/*2) x(I,_^)p7 

C2] P-*-(/*2)xP 

L3] J,£)p ( («?e) + .x 

7 



7,VP.Vi/C0]7 

7 MEliU RENU GOTOiXilUMiNUMxII 
Cl] ’ ’ 

C2] u-.'TYPE THE UUM3ER OF THE PROGRAMiS) YOU >^’ANTr 
l3] 

[4] .VP. 7^ 

C5] xuu piJu:-i^,\H 
L3] I 1^0 

C7] LQi^UUUM<II^II^^) /L 

C3] il:iIP * J J ] € ( p MEi'IU ) Cl]t* *0123456739* * ’ THEN '*L2' ELSE 

C9] L2:tIF NUMLin THEN GOTO 
[10] -^£,0 
7 
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■U-lAKEBei\Z'[^ 

V /iAKEBE 

[ 1 ] r-^£x«j(£.r,^) pP 

L2J (£./,«/ )p_d+£=0 

y 



v,-// 5 A:£:^’PCL']y 
V .‘I A KEEP 

Ll] P^(/.£)p£-^(£.J,£)pJ^O 

[2] LI :-( /O 

[3] PLJ; ]^({;+.xS^SCJ; ])i(C+.x5?(;^£C ;/; ]) 

[4] -LI ,p£[;J;]-(;x(5?(£.^)pPCJ;])t(£,£)pS+5 = 0 
l 5 ] 

7 



V.'!AKEPGin']V 
V AAKEPG 

[1] P-«?+/[3] r-Fx(x,J.£)p£ 

l2] £-£t«?(£, J,£)pP+p=0 

7 



7.V/^2’J<"/CU]7 
V Z^i-lATItl IJ;IiT 

Ll] fl r//J5 PROGRAM IS USEFUL FOR MATRIX INPUT . THE SYNTAX IS : 
L2] fl X^ MATIN 5 8 

[3] X IS THE MATRIX TO BE FORMED AND HAS DIMENSION (5x3). 

L4] fl THE PROGRAM PROMPTS FOR THE REST. 

L5] 2-iJ-O 

L6] L1:-(J./L1]<I-I+1)/L4 

L7] L2 ; 'INPUT ROM ' ,rl 
[a] -( J.Z[2] = pr-a)/L3 

[9] 'INPUT ERROR'. EACH ROW HAS ’ , ( ▼IJ’C 2 1 ) , ' ELEMENTS' 

[]0] -L2 

Cll] L3:2-2,r 
L12] -L] 

[13] L4:2-J«Zp2 
7 



V;’';£:i 7 i/La ]7 

7 .’-/g/ZP GOTO iXNUMxNUMill 

cn • ' 

[2] L'.'TYPE THE NUMBER OF THE PROGRAMiS) YOU WANT.' 

[3] ' ’ 

l4] menu 

L 5 ] X N JM-'-p N U M *- , J 

L 6 ] IJ -0 

L 7 ] LO :-*{XNUM<II*■II+^ ) /L 

L8] LltiiP * NUMjll ] e ( p MENU )[l]t* *0123456739’** THEN ’-L2’ ELSE *-L0’ 
l9] L2:tIF NUMilll THEN GOTO 
[10] -LO 
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V)v.4J:r5CJ3v 

V STAGE >JAITS C*0//' 1 ; J ;,7 ; J ;;-/1 ;A'2 ; >v'1 ; J 
Li] c/«^( p/il^x\CO//) [2] 

l 2] w1«^(</p1 ) , 

[3j .V2-/41.C1] </pO 

L4] ,0 ) 4-,V2 + (.V1 -;^2) t2 

L5] 5r^^;£;^o 

C 6 ] £_•*■ (0,«/)p\0 

l7] L\i*{{QSTAGE)<I^I^^) /Q 

[3] ^L1,pV^V,[1]+/(52’/^(;ECI-1 ].0) + ( {STAGEil^.J)^-^^ ) 

V 



7^c/[G]V 

7 

LI] Z^DG,jl*2 
7 



V£LU]7 

7 

Ll] y-^(PgQ) T( ptfJ »lQQ)xp-t-p=Q 

V 



7j£5r^CU]7 

7 

LI] >^^CDIAGtSt^ 

7 

7PlLU]7 
7 Z-^Pl 

[1] Z^5)(^.Jx£)p£x<9(£,J,^)p^ 

7 



7_£CG]7 
7 2-£ 

LI] 2-^PC;+/L 1 ] (^x£)*2 
7 



7PCU]7 
7 Z-^_rt 

L 1] Z'^+ZL 1 ] Px^x£*2 
V 



7£CU]7 

7 

l 1] ^( r//£T4 r1 00) xgf g=o 

7 
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